dc-Response of a Dissipative Driven Mesoscopic Ring. 
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The behavior of the dc-component of the current along a quantum loop of tight-binding electrons 
threaded by a magnetic flux that varies linearly in time <1>m it) — $f is investigated. We analize the 
electron transport in different kinds of one-dimensional structures bended into a ring geometry: a 
clean one-dimensional metal, a chain with a two-band structure and a disordered chain. Inelastic 
scattering events are introduced through the coupling to a particle reservoir. We use a theoretical 
treatment based in Baym-Kadanoff-Keldysh non-equilibrium Green functions, which allows us to 
solve the problem exactly. 
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I. INTRODUCTION 

The impressive development of nanoscience places the 
detailed understanding of quantum transport in meso- 
scopic systems among the main challenges of condensed 
matter physics. A rich variety of devices and struc- 
tures, including simple metallic wires, as well as com- 
plex molecules, where electrons are driven by some ex- 
ternal force, are the subject of experimental and theo- 
retical investigation. The driving field can be established 
by attached leads at different chemical potentials, a mag- 
netic flux when the system is bended into a ring or time- 
dependent fields. Transport properties, like the conduc- 
tance of the system, depend strongly on its microscopic 
details but they may also depend on the underlying driv- 
ing mechanism. 

A simple quasi one-dimensional annular system 
threaded by a magnetic flux is one of the paradigmatic 
devices to discuss the fundamentals of quantum effects 
dominating the electron transport. Several outstanding 
experiments Q, |^ |^ and theoretical works ^3, ^ have 
been devoted to the investigation of phenomena related 
to the Aharanov-Bohm effect in rings enclosing a static 
magnetic field. There is, however, relatively less litera- 
ture related to the case where the field changes in time. 
Among the latter category of problems, a very interest- 
ing example corresponds to that of a magnetic flux with 
a linear dependence on time: <i>Af(i) = ^t, which induces 
a constant electric field along the ring. 

This problem was discussed by Biittiker, Imry and 
Landauer in the early times of the theory of quantum 
transport 0. Their motivation was to search alterna- 
tive schemes to calculate the conductance of a system 
to that proposed by Landauer Q, where the sample is 
placed between two reservoirs at different chemical po- 
tentials. As first discussed in that seminal work, such 
device isolated from the external world, is not appro- 
priate to observe a dc-current since the electrons inside 
move coherently displaying Bloch oscillations and pro- 
ducing a pure ac-response. In a subsequent work Lan- 
dauer and Biittiker Q showed that when a dissipative 



mechanism is added to that system a dc-current is es- 
tablished. Biittiker also showed that a concrete element 
to introduce inelastic scattering events and dissipation 
is a lead connecting the loop to a particle reservoir j^. 
Soon later, Lenstra and van Haeringen suggested that a 
dc-current can be generated in that device by elastic scat- 
tering introduced by weak disorder • The possibility 
of "resistive" behavior originated in pure elastic scatter- 
ing processes generated a series of interesting discussions 
and criticisms [H El US El El El [13 • Most of the 
ideas of these works are based on adiabatic descriptions 
where the energy levels of the ring define minibands in 
a parametric representation as periodic functions of the 
flux. Within that adiabatic framework, scattering pro- 
cesses form small gaps between the so defined minibands 
while the time dependence of the field gives rise to Zener 
tunneling across them. In Refs. [isl ll^ interesting ar- 
guments emphasizing the concept of localization in the 
energy space as a consequence of disorder have been pro- 
posed against the possibility of dc-response without dis- 
sipation. 

The discussion of the role of dissipative effects in 
the driven ring was introduced on the basis of a phe- 
nomenological equation of motion that describes the re- 
laxation via inelastic scattering processes (ISP) of the 
time-dependent occupation of a miniband Q . The main 
argument was that the dc-response should vanish in the 
limit of vanishing and strong relaxation, while it should 
peak at some intermediate regime. For small electric 
fields, the behavior of the current is also found to follow 
a linear, ie. Ohmic-like, dependence as a function of the 
induced electromotive force (emf ) . These notable predic- 
tions have not been examined in more detail during many 
years. Quite recently, the effect of inelastic scattering in 
the dc-behavior of that system has been studied Ell ■ In 
that work a pure "clean" system is considered, (ie with- 
out any kind of elastic scattering). Non-Ohmic behavior 
has been found in the dc-current vs emf characteristic 
curve within the limit of small ISP, while the dc-current 
is found to decrease continuously as the strength of the 
ISP increases. The fact that no tendency towards a van- 
ishing component of the dc-current is observed as the 



dissipation tends to be suppressed is rather surprising, 
since the limit of pure Bloch oscillations in the isolated 
system seems not to be recovered. 

In this work we consider a ring threaded by a flux with 
a linear time dependence and we analyze in detail the 
role of dissipation in that system. We study the clean 
one-band system as well as the effect of two different 
kinds of elastic scattering processes: a periodic poten- 
tial with a two-sublattice profile and a random potential. 
The effect of inelastic scattering is modeled by means of 
the coupling to a particle reservoir through an external 
lead. We use a recently proposed theoretical treatment 
based on Baym-Kadanoff-Keldysh Green functions 
which enables a full out-of-equilibrium description of the 
time-dependent problem, without introducing any kind 
of adiabatic assumptions or approximations. As far as no 
many-body interactions are considered, that treatment 
leads to the exact solution of this problem. Thus, one of 
the aims of this work is to analyze the case of the dis- 
ordered driven ring in the limit of weak ISP, which has 
focused several efforts but has never been tackled with 
an exact method. The formalism based on Green func- 
tions is particularly appealing to deal with the coupling 
to external leads and reservoirs. These elements can be 
introduced by means of physical models and they can be 
concisely represented by self-energies, avoiding further 
assumptions about the boundary conditions. We con- 
sider two different electronic models for the lead: one is 
described by a constant density of states and a very wide 
band, while the other is represented by a semi-infinite 
tight-binding chain. We show that the dc-response does 
not depend qualitatively on the latter details. In all the 
cases, our results indicate that the dc-current tends to a 
vanishing value as the coupling to the reservoir goes to 
zero, displays a maximum and then decreases as that cou- 
pling becomes stronger. For clean systems, the position 
of the maximum shifts towards very low ISP strengths as 
the emf decreases and as the length of the ring increases. 

The system considered here, actually lies in the cate- 
gory of the so called ratchet problems, since the induced 
electric field by itself is not able to produce a dc-response, 
needing the aid of some additional rectification mecha- 
nism. This motivates the analysis of the transport prop- 
erties of this system in the framework of recent symmetry 
arguments proposed to examine rectification processes in 
ratchet systems in the presence of time-periodic fields 
|20l I21I l2a | . Such arguments provide further support to 
the idea that dissipation is an essential ingredient for this 
system to have a dc-current irrespectively the particular 
nature of the elastic scattering along its circumference. 

The paper is organized as follows. In section II we 
present the model and the theoretical approach to ana- 
lyze the behavior of the current. In section III we present 
the results. Section IV is devoted to summary and con- 
clusions. 




FIG. 1: (Color on line) Scheme of the setup. The system 
indicated in the red box defines the reservoir with chemical 
potential fj,a ■ 

II. THEORETICAL TREATMENT 
A. The model 

We consider the device sketched in Fig. It consists 
in a ring threaded by a magnetic flux with a dependence 
on time of the form ^M{t) — ^t, in contact to a particle 
reservoir with a chemical potential /i^ through an exter- 
nal lead. The full system is described by the following 
Hamiltonian: 

H = Hring{t) + Ha + i?la, (1) 

where the first term, representing the ring, depends ex- 
plicitly on time due to the presence of the time-dependent 
flux. We consider noninteracting spinless electrons de- 
scribed by a tight binding model with N sites and lattice 
constant a = L/N. We also consider the possibility of 
elastic scattering in the system, described by a local en- 
ergy with a profile e/ with I = 1, . . . , iV. The Hamiltonian 
is 

N N 

1=1 1=1 

.(2) 

with the periodic condition iV -I- 1 — 1. The time- 
dependent phase (jft attached to each link, with (f> = 
^/{^qN), being $0 = hc/e, accounts for the presence 
of the external magnetic flux. 

The term Ha describes the reservoir. We consider 
two models of noninteracting electrons with a band- 
width W and a chemical potential fia for this system: 
(i) a wide-band model, defined by a constant density 
of states and a very large W, and (ii) a semi-infinite 
tight binding chain with hopping amplitude IF/ 2, which 
corresponds to a semicircular density of states pa{^) = 
4e(|w| - W)y/W'^ - uo"^ /W"^ . 

The last term of H , 

Hia = -Ta{c\ca + cJ^Ci), (3) 
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represents the connection between the ring and the reser- 



B. Dynamical equations and symmetry properties 

In Refs. [20I I21I an interesting connection has 
been suggested between the rectification properties of the 
system and the underlying symmetries of the equations 
of motion in several ratchet problems. The main idea 
is that the symmetries of the equation of motion that 
change the sign of the velocity would lead to a vanish- 
ing dc-component of the current. We now turn to follow 
the lines suggested in those works to carry out a similar 
symmetry analysis in our problem. 

Let us first consider the simpler case of the ring isolated 
from the reservoir, the problem is described by Hring 
alone. For sake of simplicity in the notation we shall 
adopt a system of units where $0 = 1, H = 1. The 
relevant equation of motion is the Schrodinger equation 
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eik{t)ipkit), 



(4) 



being eik{t) = -Th{e-"^Hk,i+i + e"^Hk,i-i). We have 
adopted the Einstein summation rule over site indexes. 
The current in the isolated ring results 

JtT-liit) - eT,Im[e-^*VrW^;+iW]- (5) 
The matrix elements of the Hamiltonian satisfy 



'Ik 



{-t)=eik{t) 



-kl 



[t). 



Hence, the inversion t ^ —t followed by the complex 
conjugation of Q leads to ipi{—t) = ipi{t) which, when 
replaced in ij^l, results in the following property of the 
current 



jisol / *\ jisol fi\ 



(7) 



Since in our reasoning we have not assumed any partic- 
ular form for the energy profile ei, the criterion of Refs. 
[20I I21I rules out the possibility of a dc-current as 
far as time-reversal symmetry is preserved in the device 
irrespectively the particular model assumed for the elas- 
tic scattering processes. Another symmetry operation 
that changes the sign of the current is spatial reflection, 
however, the latter is not a symmetry of the Schrodinger 
equation 10}. 

In the more general case, when inelastic scattering pro- 
cesses are considered, the relevant equation of motion is 
the Dyson equation for the Green function. In the frame- 
work of Baym-Kadanoff-Keldysh formalism the latter has 
a matrix form and one must work with two independent 
Green functions: the retarded Green function, 

G5(t,O=-*0(i-O([c.W,c](i')]>, (8) 



and the lesser Green function 



(9) 



The latter determines the mean values of the observables. 
In particular, the current through a bond {1,1 + 1) is 
written as 



(10) 



The Dyson equation for the matricial Green function 
leads to the equations of motion for the retarded and 
lesser components 

-^±G<it,t')-G<{t,t')sk,{t') = 
dt"[Ggit,t")^<it",t') 



+G<{tX)^tS",t' 



(11) 



while the evolution for the lesser component at equal 
times is given by 

dG<(t,t) 

' dt " ^^^""^ + ^■^kmt^ {t, t) 

= j dt''[^git,nG<it'',t)+^<it,nGt^it'',t) 

-Gf,(t,t")S,< - G< (t,t")S^,r,t) ] , (12) 



(6) being G^^{t,t') = [Gfjt',t)]* and S^^.(<,t') = 



In the Hamiltonian limit, ^^jit,t') = (i,t') = 0, 
in the absence of many-body interactions. Making use of 
the property of the matrix elements of the Hamiltonian 
and of the symmetry property of the Green function 
Gf/t,t') = -[G,< (t',t)]*, the property for the time- 
inversion operation on the current is recovered. 

It is easy to verify that any non-vanishing self-energy 
correction that fails to satisfy the properties S*^,(i,t') — 
^Ut,t') = 0, Sffc(i,f) = s,k{t)S{t-t'), with s„(t) = 
— breaks time- reversal symmetry in the 

equations of motion. In particular, any self-energy rep- 
resented by kernels of the form 



^git,t') = -iQ{t-t') 



dcu 
2^ 



{-2lm[j:g{oj)]}e 



j(t-t') 



^fkM 



duj 
2^ 



/(c.){-2Im[Ef,(c.)]}e— (*-*'), (13) 



with Im[I]^(aj)] 7^ 0, being /{uj) the Fermi function, will 
break time-reversal symmetry in the equations of motion 
for the Green functions and will break the property ^ 
for the time-dependent current. 

In our problem, the effect of the exchange of parti- 
cles and energy between the mesoscopic system and the 
reservoir can be exactly written in terms of a self-energy 
correction at the site / = 1 of the tight-binding ring 
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[l^ 0, • The retarded and lesser components of this 
self-energy are 

= -ze(i-t')|T„P 

= (J^a^.iSfli-i') (14) 
= ^\Ta\'J ^/(^)Pa(^)e-'"(*-*')5.,i4,i, 

where Pa{uj) is the density of states of the reservoir 
and the Fermi function is /(lu) = i/(e'3("-'^°) + 1). 
In our calculations we consider zero temperature, i.e. 
f{uj) = Q{uj — Ha)- The wide-band model leads to a con- 
stant imaginary retarded self-energy, (w) = ia, and 
S^(tj) = ifa{uj)a, being a — \Ta\'^TT/W. The dissipative 
nature of the coupling to the reservoir manifests itself in 
the fact that 'Si'{u!) has a finite imaginary part. As dis- 
cussed above, this breaks time-inversion symmetry in the 
dynamical equation, removing the inversion of the cur- 
rent under this symmetry operation. This satisfies the 
criterion of Refs. [20l I21I l2a |. namely, a non- vanishing 
net current is possible when symmetries of the equation 
of motion leading to an inversion of the time-dependent 
current are broken. 



C. Evaluation of the Green functions and the 
dc-current. 

The formalism leading to the evaluation of the cur- 
rent along the ring has been presented in Ref. ^^I- In 
this subsection we summarize the main equations and we 
defer the reader to this work and references therein for 
further details. 

Following Ref. 19] , it is convenient to perform a gauge 
transformation in the fcrmionic operators of the Hamil- 
tonian 

c„ = exp[m0i]c„, (15) 

according to which the Green function for the positions 
m, n on the ring transforms as 

G^Jt, t') = eM^^{mt - nt')]GZn(.t, t'). (16) 
Defining the Fourier transform 

GtnitX) - / ^G^,„(t,^)e— (*-*'), (17) 
the resulting equation for the retarded Green function is 
Gf,,„(t,w) - G^^„(u;)-G,%(i,c. + <i>)Tjvi(OG°„(o.) 



with 

G^,i(t, u) -f Gl^j,{t, u + ^)Tni ml.ico) 
+Gf„^,{t,Lu- $)riA,(OG°.,i(c.) = G?„,i(c.) 
G^^^it, u:) + gI^mH, ^ + '^)TNi{t)G%{io) 
+G^,^,{t,cu- <^)T,N{t)G%^^{Lj) = Gl,M, (19) 

where TiN{t) = [Tmit)]* = Tue^'^K For each time t, the 
solution of the above set of linear equations provides the 
complete exact solution of the problem. 

The equilibrium Green function G^„„(ti') corresponds 
to the problem of an open chain in contact to the reser- 
voir. It is obtained from the solution of the Dyson equa- 
tion 

Gl^ni^) = g°.,„(^) + G^^i(c.)Sf (c.)g?,„(c.), (20) 
where 

^ 1 

9°m ni^) KAl F (21) 

i^— 1 

is the Green function of the chain isolated from the reser- 
voir. It can be expressed in terms of the eigenvalues E^, 
and eigenvectors — J^i ^rlO of the Hamiltonian 

N-l N 

Ho = -n J2 C^h+i + cJ+iQ) + J2{Vi + ei)clci, (22) 
1=1 1=1 

where Vi — (t>l is the scalar potential due to the induced 
electric field. 

The set (|18ll and (|19|) describes the process of the clos- 
ing of the biased chain by means of an effective time- 
dependent hopping Tijv(i) that pumps electrons with a 
frequency $ through the bond (1-/V). This set involves, 
in principle, an infinite number of equations. In the nu- 
merical procedure, upper and lower energy cutoffs Ai, A2 
are chosen such that Ai < < A2, '^v. 

Starting from its definition pUI) , the current along the 
ring can be written 



JLi+i{t) = 




xG,+i,i(i,c. + 0)[G;,i(t,c. + </>)]*], (23) 

where we have used G^i{t,t') = [Gf^^{t' ,t)]* . Therefore, 
the evaluation of the time-dependent current is reduced 

to the evaluation of Gi i{t,uo) obtained from ifTiUl . This 
quantity displays an oscillatory behavior as a function of 
t, with the period tb = 27r/$ of the Bloch oscillations. 
The dc-component does not depend on the bond {1,1 + 1) 
chosen for the calculation and it is defined as 

Jdc^— I dtJi,i+i{t). (24) 

TB Jo 
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III. RESULTS 

This section is devoted to analyze the behavior of the 
de-current H24|l as a function of the induced emf $, the 
strength of ISP and the chemical potential The 
strength of ISP is related to the degree of coupling be- 
tween the ring and the reservoir. In the case of the wide 
band model, the parameter cr sets that measure. In the 
model with a semicircular density of states, that mea- 
sure is given hy T^/W. In our calculations, we fixed 
W = ATfi and changed Tq,. All energies will be expressed 
in units of the hopping parameter Th- We shall anal- 
ize three different energy profiles ei for the tight-binding 
model that define a clean one-band chain, a chain with a 
two-sublattice structure and a random potential. 

A. The clean one-band ring 

This case corresponds to e; = 0, Z = 1, . . . , in ||2Jl. In 
the limit of vanishing dissipation (the ring isolated from 
the reservoir) the system is described by the Hamiltonian 
Hring{t) alouc, the problem has time-reversal symmetry 
and a vanishing dc-current is expected. In fact, for ei = 0, 
Hring{t) Can bc casily solved by performing a Fourier 
transform to the fc-space. The retarded Green function 
can be calculated analytically, resulting G^ „(<,t') = 

— — t') ^ji, exp{jfc(m — n)} exp{— j (isefc(s)}, with 
e-k{s) = — 2T/iC0s(fc -I- 0s), being k = 2m:/L,n — 
0, . . . ,N. The current can also be obtained analytically, 
resulting 

J{t) = Th sin(fc + (j}t) . (25) 

k 

where the summation extends over the set of occupied 
fc-states. 

The behavior of the dc-current in the coupled system 
at a fixed chemical potential /Iq for different strengths of 
ISP and $ is analyzed in Fig. |21 A wide band model was 
assumed for the reservoir in these calculations. The dc- 
current displays a maximum which shifts to lower values 
of (7 as $ decreases. At low cr, the imaginary part of 
„(w) evolve towards a sequence of delta functions and 

the Green functions „(t, uj) also develop a structure of 
peaks that get sharper. From the practical point of view, 
this turns harder the evaluation of the integral in ui (|23|l . 
preventing us from exploring the range of cr/T/j < 0.05. 
However, we think that the range considered is enough to 
infer the trend towards cr — s- 0. The characteristic curves 
Jdc vs $ are shown in the inset for particular values of 
a within the regimes with low, intermediate and strong 
coupling to the reservoir. In the latter limit, the dc- 
current remains linear within a wide range of $ while for 
weak coupling, the departure from linear response takes 
place at small $. 

To study the origin of this peculiar behavior we must 
analyze the structure of the equations (|18|l and H19|) . At 
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FIG. 2: (Color online) dc-current as a function of cr in a ring 
with N = 20 sites for different values of the induced elec- 
tric field -J- = 0.01,0.02,0.04,0.06,0.08,0.1,0.2,0.3,0.5 (up- 
per panel) and with N = 200 for $ = 0.2, 0.3, 0.4, 0.5 (lower 
panel). The chemical potential of the reservoir is ~ — 1- 
The plots corresponding to the lowest and the highest <E> are 
drawn in white and black symbols, respectively, while grey 
symbols correspond to intermediate values. Inset: Jdc vs $ 
characteristic curve for different strengths of ISP correspond- 
ing to cr = 0.05, 0.2, 0.8 (red circles, green squares and blue 
diamonds) . 



each time t, the Green function G^^{t,iu) contains a 
combination of a large number of w-components (sep- 
arated in Auj — <&) of the Green function (|20() . In 
the weak coupling limit, G^,j„(w) is sizable only within 
a neighborhood of the N frequencies E^. Therefore, 
a resonant combination of components at different fre- 
quencies is achieved when Auj ~ AE^, the latter being 
the mean energy separation between two eigenenergies, 
AE^ = {E^+i - E^) - 4Th/N. Whether such interfer- 
ence will be constructive or destructive and give rise to 
a large or a small component of the dc-current is a ques- 
tion without obvious answer, particularly in the case of 
small $, where the number of coupled frequencies is very 
large. On general grounds, one expects that such an ef- 
fect would strongly depend on the underlying symmetries 
of the model. Related discussions have been recently pre- 
sented in problems of pumping introduced by microwave 
fields [2^|27n2SiJ where the Floquet representation of the 
wave function is used, leading to a structure of the so- 
lution containing a mixing of frequencies that differ in 
the frequency of the pumping field. In our problem, it 
becomes apparent that such interference is constructive 
regarding the magnitude of Jdc and for this reason, in 
the weak coupling limit, the current grows linearly as a 
function of reaching a maximum at $ ~ AE^. The 
decrease for $ >> AEy can be understood by noting 
that the number of resonances is approximately 4Th/^, 
becoming smaller as $ increases. The maximum of Jdc 
as function of a can be also interpreted in terms of com- 



FIG. 3; (Color online) dc-current as a function of the strength 
of ISP for a model with a semicircular density of states 
with bandwidth W = 4. Different plots correspond to 
$ = 0.1,0.2,0.3,0.4,0.5. The plots in white and black sym- 
bols correspond to the lowest and the highest "1>, while the 
grey ones correspond to the intermediate values. The chem- 
ical potential of the reservoir is = —1. Inset: Detail of 
the behavior of Jdc for low ISP. Red circles and blue squares 
correspond to $ = 0.4, 0.5, respectively. The behavior is con- 
sistent with Jdc tx Ta when ^ 0. 



bination of components of the Green function from dif- 
ferent frequencies. In fact, the effect of increasing ISP 
is to spread the spectral weight of the peaks of ni^)- 
Thus, for strong couphng the interference involves a large 
number of frequencies but with a low (approximately con- 
stant) weight. For lower ISP, the amplitudes can be large 
for frequencies close to Ei, but tend to be vanishingly 
small in between. The result is that there is a maximum 
in Jdc at some strength of ISP that seems to scale as 
oc at least, within the regime $ > AE^. On the side 
of strong coupling, where the resonant effects are highly 
smoothed, the resulting dc-current is relatively smaller 
(consistent with the idea that resistance increases with 
ISP) and behaves linearly within a larger range of $. 

The effect of the length of the chain is also analyzed 
in Fig[21 The first issue to note is that in the N — 200 
sites ring (see lower panel of Fig. 1, the dc-current for 
$ = 0.2 remains growing for the smallest a considered. 
Instead, for = 20 and the same parameters, the cur- 
rent decreases down from its maximum (cf upper panel 
of Fig. 121) . This behavior is consistent with the picture 
of constructive resonances: The condition (f> > AE^ is 
achieved in this case at lower $, while the larger number 
of peaks increases the probability of resonances and the 
current remains large, within a larger range of <&. For 
larger fields and for larger tr, the qualitative and quanti- 
tative behavior is essentially the same as that observed 
in the smaller ring. 

Regarding the effect of the model for the reservoir, it 
becomes clear from the results of Fig. O that it does 



FIG. 4: (Color online) dc-current as a function of the chemical 
potential of the reservoir (represented by a wide band model) 
for two values of the emf <1> — 0.1, 0.5 (solid and dashed lines) 
and two strengths of ISP a — 0.05, 0.85 (thick black and thin 
red lines), respectively for a ring with — 200 sites. 



not play any relevant qualitative role since we can iden- 
tify in these plots the same features observed in Fig. |2 
In Ref. |2^ an interesting connection have been found 
between the behavior of Jdc as a function of when 
—> and the underlying symmetries of the Hamil- 
tonian. The plots of Fig. corresponding to the two 
highest emfs ($ — 0.4, 0.5), where the decrease to a van- 
ishing can be cleanly captured wihin the shown range 
of T^, suggest a Jdc oc behavior as T„ —^ 0. This 
is even more clear in the plot of Jdc/T^ vs shown in 
the inset. This dependence can be understood by noting 
that when AE,, + ^ » T^/W, the retarded Green func- 
tions entering the evaluation of J/_i+i(t) (see Eq. E^ . 
can be replaced by the ones for the uncoupled ring (cor- 
responding to = 0). Hence, for small enough T^, the 
time-dependent current and also the dc-component Jdc 
become linear in T^. This is in perfect agreement with 
the analysis done in Ref. |26| for a linear device pumped 
with a laser in configurations with broken time-reversal 
symmetry. 

The behavior of Jdc as a function of the chemical po- 
tential is shown in Fig^ It is remarkable that a semi- 
circular shape can be identified in the envelope of these 
functions. As mentioned before, the density of states of 
a semi- infinite tight-binding chain of bandwidth W , is 
proportional to a semicircle of radius W/2. This moti- 
vates the comparison of the conductance of our device 
with the conductance of a quantum dot which has two 
semi-infinite tight-binding chains with hopping elements 
Tfi and chemical potential jja attached to its left and to 
its right, through a hybridization amplitude Weff- The 
latter corresponds to the usual configuration employed in 
a Landauer-like calculation of the conductance through a 
dot connected to a left and a right lead. Assuming that a 
small bias V is applied between the left and the right, the 
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FIG. 5: (Color online) F (proportional to the effective density 
of states peff{fJ,a) obtained by assuming a model for the con- 
ductance based on a dot coupled to two semi-infinite chains 
for the case of a ring with TV = 20 sites and a wide band 
model for the reservoir. Thick and thin solid lines correspond 
to a = 0.05 and a — 0.85, respectively. Black and red lines 
correspond to $ = 0.01 and <& — 0.1, respectively. The plot 
in blue dashed lines corresponds to "1> = 0.1 and a model of 
reservoir with a semicircular density of states with W = 4 
and = 0.1. 



all the features in the weak coupling regime observed 
within the wide band model for the reservoir are also 
obtained with a reservoir with a semicircular density of 
states (cf plot in dashed hnes of Fig. |3J). In the strong 
coupling limit, Weff ~ 3 is approximately the same for 
both fields and a fine structure of N wide peaks is dis- 
tinguished in peff{fia)- The effective parameters F and 
Weff do not have any straightforward significance in the 
context of our original model and the parallel between the 
driven ring and the dot connected to semi-infinite chains 
is purely heuristic. However, it is remarkable the similar 
behavior of the conductance of the two devices. It sug- 
gests that although the chain that forms the ring is made 
up by a discrete chain with a finite number of sites, the 
dressing due to the mixing of a large number of frequen- 
cies tend to produce a structure such that any site along 
the ring would feel as if placed between two semi-infinite 
leads. In spite of this resemblance, the quantitative prop- 
erties of the conductance are different from the one which 
would be measured in a Landauer device with the chain 
forming the ring stretched and placed between two reser- 
voirs, as already discussed in [S^l. In fact, neither the 
effective density of states at the site between the two ef- 
fective leads nor the hopping element Wef / correspond to 
the configuration of a non-interacting site coupled to two 
leads through as in the original chain. 



conductance of such a system is given by a Landauer-like 
formula [2^ 

Geff = = epo(/^a)r, (26) 

where poif^a) = ^4T,^ — is the density of states 

of the semi-infinite chain and F — Wgj^pe//(MQ)/47r, be- 
ing peff{fJ.a) the density of states of the central system 
dressed by the contact with the two chains. 

In order to define the conductance for the device of 
Fig. n let us note that the induced electric field is E — 
^/(Lc). Therefore, if we make the standart assumption 
that for small the current can be written as Jdc = f^cE, 
being CTc the conductivity, the conductance of the chain 
is defined from Jdc = GEL. Therefore, in units where 
(c = 1), the conductance for small enough bias is 



If we adopt the expression H2t)|l as a phenomenological 
model for the conductance of our system, we obtain for 
F the result shown in Fig. [S] In the limit of weak cou- 
pling, this behavior is consistent with an effective den- 
sity of states approximately constant along the band- 
width. Taking into account its sum rule, it should be 
PeffifJ-a) ~ <d{\fic\-4:Th)TT/'iTh. The large average value 
of F indicates a large effective hopping from the dot to 
the semi-infinite chains {weff ~ 12,8) for $ = 0.01,0.1, 
respectively. The details of the model used for the reser- 
voir do not influence the final behavior of F. In fact. 



B. The ring with a two-band structure 

As mentioned in the introduction, one of our motiva- 
tions is to analyze the behavior of the dc-current in a 
disordered system in order to explore the possibility of 
'resistive' behavior in the limit of vanishing dissipation. 
We found it instructive to study first an intermediate sit- 
uation where the potential profile is e2n = Ce, ^2n-i — 
eo,n = 1, . . . , N/2. This structure defines two bands sep- 
arated by an energy gap which are the basic ingredient 
to discuss the effect of Zener tunneling. 

The dc-component of the current as a function of a 
is shown in Fig. for a profile with = 0.2, Co = 
and chemical potentials of the reservoir within the lower 
band (upper panel) and within the gap (lower panel). 
In the first behavior similar to that found in the 

clean one-band model (cf Figs I2I3|I is observed. Namely, 
a maximum of the current that shifts to lower a as $ 
decreases. The shift depends, however, slower with $ in 
the present case: taking as a reference the plot corre- 
sponding to $ = 0.2, we see that the maximum in the 
two-band case takes place at cr '--^ 0.2 while in the clean 
ring Jdc remained growing for the lowest a used in our 
calculations. The comparison of the plots of the upper 
panel of Fig. [7| with those of the right panel of Fig. O 
shows that the magnitude of the maxima in the weak 
dissipation limit is smaller in the two-band ring than in 
the case of the clean ring with the same number of sites 
in contact to a reservoir with the same characteristics 
and with the same chemical potential. The behavior of 
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FIG. 6: dc-current as a function of a for two different chem- 
ical potentials in a ring of = 200 sites with a two- 
band structure defined by Co = 0, te = 0.2 and a reser- 
voir with a constant density of states. Different plots cor- 
respond to $ = 0.05,0.1,0.2,0.3,0.4,0.5 (upper panel) and 
$ = 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8 (lower panel) . 



Jdc as a function of the chemical potential shown in Fig. 
[3 shows that this is also the case for other values of the 
chemical potential and dissipation strengths. The smaller 
current for chemical potentials within the lower conduc- 
tion band, indicates that the presence of an energy gap 
in the structure of energy levels of the chain contributes 
to a less efficient mixing of weights of the Green function 
at different frequencies relative to the corresponding one 
in the clean one-band case. 

The effect of ISP on the magnitude of the current 
through the gap is shown in the lower panel of Fig. 
Within the weak coupling regime, the current is vanish- 
ingly small even for fields significantly larger than the en- 
ergy gap (~ |eo — Eel). This indicates that Zener tunneling 
alone (ie without dissipation) is not enough to generate 
a dc-response. 

The behavior of Jdc as a function of the chemical po- 
tential is shown in Fig [T] An important asymmetry is 
observed between the upper and lower bands, which is 
more pronounced for low fields. This is an indication 
of the relevance of the interband processes generated by 
the coupling of the electrons with the field. This discour- 
ages us from following the steps of the previous section 
in trying to make a parallel with a device based in two 
semi-infinite leads with the same band structure of the 
ring, since the latter would have two symmetric bands. 



C. The ring with a random potential 

This final subsection is devoted to analyze a disor- 
dered ring, with a random potential e; = twli, being 
^1 < 7/ < 1, a random number. The dc-current as a 
function of ISP strength for a fixed chemical potential 
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FIG. 7: (Color online) dc-current as a function the chemi- 
cal potential for a ring with = 200 sites and a two-band 
structure defined by to = 0, £e = 0.2 and a reservoir with 
a constant density of states. The upper and lower panels 
correspond to weak {a — 0.05) and strong (cr = 0.85) cou- 
pling regimes. Different plots correspond to different fields 
$ = 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8. The ones corresponding 
to the lowest and the highest fields are plotted in thick black 
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FIG. 8: Color online) dc-current as a function of the strength 
of ISP for a ring of A'' = 200 sites with a random potential with 
amplitude = 0.2 and a reservoir with a constant density of 
states. Different panels correspond to different realizations of 
the random potential. Black circles, red squares, green dia- 
monds and blue triangles correspond to $ = 0.2, 0.3, 0.4, 0.5, 
respectively. 
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FIG. 9: (Color online) dc-current as a function of the chemical 
potential for a ring of N = 200 sites with a random potential 
with amplitude e^, = 0.2 and a reservoir with a constant den- 
sity of states. The upper and lower panels correspond to weak 
(cr = 0.05) and strong (a — 0.85) coupling regimes. Different 
plots correspond to different fields $ = 0.1,0.2,0.3,0.4,0.5. 
The ones corresponding to the lowest and the highest fields 
are plotted in thick black and violet lines, respectively. 

is shown in Fig. |S1 for four different realizations of the 
random potential. In all the cases a similar behavior is 
found: the current exhibits a wide and mild maximum 
at intermediate coupling and tend to zero in the limit 
of vanishing ISP. The magnitude of J^c at the maxima 
is smaller than in the two-band ring and significantly 
smaller than in the clean one-band system. As mentioned 
in the previous subsection the small current in the limit 
of weak coupling to the reservoir can be interpreted by 
noting that the formation of energy gaps in the struc- 
ture of the energy levels of the chain tend to break the 
resonant behavior presented in the clean chain for small 
fields in the weak coupling regime. For small enough cr, 
the trend is Jdc oc cr, as in the clean limit. 

As a function of the chemical potential, the dc-current 
displays fluctuations that become stronger in the strong 
coupling limit (see Fig. O . The fluctuations bear a close 
resemblance with those observed in the behavior of the 
conductance of disordered quantum wires [sj. A sys- 
tematic analysis of the probability distribution of con- 
ductance G defined in (|27|) over several disorder realiza- 
tions as well as the behavior of G as a function of the 
system length is left for the future. 

IV. SUMMARY AND CONCLUSIONS 

We have studied the transport properties of a dissi- 
pative tight-binding ring driven by means of a magnetic 



flux with a linear dependence in time. Dissipation is in- 
troduced by coupling the mesoscopic ring to an external 
macroscopic system that plays the role of a reservoir of 
particles and energy. We have extensively analyzed the 
behavior of the dc-component of the current along the 
ring as a function of the strength of dissipation and of 
the chemical potential of the reservoir for clean chains, 
chains with a two-band structure and an energy gap and 
disordered chains. 

We have analyzed the influence of the specific model 
assumed for the reservoir and we conclude that this does 
not play any relevant role. In the weak coupling regime, 
the structure of the energy levels of the Hamiltonian de- 
scribing the ring, seems to play a relevant role. Such 
structure determines the way in which contributions from 
different frequencies of the Green function couple through 
the pumping term in Eqs. ifTHl) . l(TI?|l . 

In all the cases we expect a vanishing dc-current in the 
Hamiltonian limit (vanishing ISP). For the case of the 
clean ring, we find that the current exhibits a maximum 
as a function of the ISP strength that seems to scale 
with For small enough $ this maximum can lye very 
close to the limit of vanishing ISP. The range of fields 
where a linear behavior of the dc-current as a function 
of $ is observed depends on the strength of ISP. It is 
wide for strongly coupled systems and very narrow for 
weakly coupled ones. Remarkably, the behavior of the 
conductance of the system as a function of the chemical 
potential can be reproduced with a device consisting in 
a quantum dot with a constant density of states and two 
semi-infinite tight-binding chains with the same hopping 
parameter as the original one, attached to its left and 
right sides through a very large hopping element. 

In the ring with a two-band structure, the behavior of 
the dc-current within each of the two conducting bands is 
similar to that observed in the one-band case. However, 
its magnitude is smaller and at fixed chemical potential 
and emf , the position of its maximum as a function of the 
ISP strength has a softer dependence with the emf. An 
important asymmetry is observed between the behavior 
of the current within the upper and lower bands, while 
inelastic scattering is essential to obtain a sizable mag- 
nitude of the current through the gap, irrespectively the 
intensity of the induced emf. 

The disordered ring can be seen as the N-band ex- 
tension of the two-band case. We have examined some 
typical realizations of the random potential and observed 
that as a function of inelastic scattering, the maximum 
in the dc-current takes place within the intermediate 
regime, becoming vanishingly small in the limit of weak 
coupling. As a function of the chemical potential, the 
current displays fiuctuations following patterns that de- 
pend on the degree of inelastic scattering. 

Our results are in agreement with the criteria based 
on symmetry arguments suggested for ratchet problems 
[20I I21I |23 | , according to which a dc-component of the 
current is possible provided that the equations of motion 
of the system are not invariant under symmetry oper- 
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ations that change the sign of the time-dependent cur- 
rent. In the case of the present problem the only pos- 
sible symmetry that may introduce such an inversion, 
is time-reversal. This symmetry is an exact one in the 
pure Hamiltonian limit where the ring is isolated from 
the reservoir but it is immediately broken when the cou- 
pling to the reservoir is considered. Furthermore, we have 
found that the current grows linearly with the parameter 
that characterizes the strength of the inelastic scatter- 
ing. This behavior has also been found in other pumped 
systems with broken time- reversal symmetry |2a |. The 
behavior of the dc-current in the case of the two-band 
and N-band (disordered) rin g is also in agreement with 
the conclusions of Refs. [H 111 El El lli which have 
been devoted to argue against the possibility of resistive 
behavior caused by Zener tunneling alone. In our case we 
were able to solve the problem exactly and to examine de- 
tails of the behavior of the current for different strengths 
of inelastic scattering. For this reason, we think that our 



results provide a robust support to the idea that dissipa- 
tion is an essential ingredient to obtain resistive behavior 
in this system. 
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